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Abstract. We investigate a calculation method for solving the Mukhanov-Sasaki equation 
in slow-roll fe-inflation based on the uniform approximation (UA) in conjunction with an 
expansion scheme for slow-roll parameters with respect to the number of e-folds about the 
so-called turning point. Earlier works on this method has so far gained some promising 
results derived from the approximating expressions for the power spectra among others, up 
to second order with respect to the Hubble and sound flow parameters, when compared to 
other semi-analytical approaches (e.g.. Green’s function and WKB methods). However, a 
closer inspection is suggestive that there is a problem when higher-order parts of the power 
spectra are considered; residual logarithmic divergences may come out that can render the 
prediction physically inconsistent. Looking at this possibility, we map out up to what order 
with respect to the mentioned parameters several physical quantities can be calculated before 
hitting a logarithmically divergent result. It turns out that the power spectra are limited 
up to second order, the tensor-to-scalar ratio up to third order, and the spectral indices and 
running converge to all orders. This indicates that the expansion scheme is incompatible 
with the working equations derived from UA for the power spectra but compatible with that 
of the spectral indices. For those quantities that involve logarithmically divergent terms in 
the higher-order parts, existing results in the literature for the convergent lower-order parts 
calculated in the equivalent fashion should be viewed with some caution; they do not rest on 
solid mathematical ground. 
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1 Introduction 

Cosmological inflation [1-8] is a rapid exponential expansion of spacetime. It stretched the 
quantum fluctuations which eventually gave birth to large-scale structures that we nowadays 
observe as galaxies and clusters of galaxies. Being able to provide an elegant solution to 
the homogeneity and flatness problems associated with hne-tuning the appropriate initial 
conditions for the Universe and explain the emergence of the seeds of structure formation, 
it serves as a good complementary patch for the standard Big Bang Cosmology. What we 
have today is a paradigm with great explanatory power involving a hybrid theory composed 
of inflation and Big Bang theory. 

Owing to its incomplete nature however, inflationary cosmology is not without hurdles. 
When making predictions or deriving expressions for physical quantities that we want to 
compare with experiments, one is faced with at least two possible problems. The first one 
deals with the possible disagreement of the prediction with experiment. The second one 
involves the self-consistency of the prediction. Focusing on the latter, one may consider for 
instance, the possible emergence of divergent quantities that need to be regularised or renor¬ 
malised (see for instance. Refs. [9-11]), contradictory ends of calculations of non-Gaussianity 
(see Refs. [12] for a review of non-Gaussianity), or some mathematical artifacts that when 
scrutinized may lead to misleading results. 

We take the last point above seriously in examining the calculation of the power spec¬ 
trum among others, carried out in Ref. [13] (see also Refs. [14-16]), using the method of the 
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uniform approximation (UA) [17, 18]. Our inquiry centres around the possible emergence 
of logarithmic divergences as mathematical artifacts of one calculation scheme for relevant 
physical quantities in inflationary cosmology. To elaborate, we start with the Mukhanov- 
Sasaki equation for scalar perturbations in the inflation model we are interested in this work 
namely, slow-roll fc-inflation [19, 20]. It states 

4' + - Us) Vk = 0, (1.1) 

where is related to the gauge-invariant cosmological perturbation TZ^ [21] as Vk = zTZk 
with and a is the scale factor, e = is the first slow-roll parameter 

defined here in terms of the Hubble parameter H, k is the wavenumber, Cg is the speed of 
sound. Us is the effective potential, and the symbols dot and prime denote differentiation with 
respect to the coordinate time t and conformal time r] respectively. The effective potential 
takes the form Ug = {l/rf‘){2 + ■ ■ ■), where the dots from the point of view of slow-roll k- 
inflation, represent quantities that are small in comparison to the addend two. This form of 
the potential is suggestive of a perturbative solution. Furthermore, owing to the structure of 
the Mukhanov-Sasaki equation pertaining to its behaviour for small- and large-wavelength 
modes, it is subject to semi-analytical approaches to solving and hence, calculation of 
relevant quantities like the power spectrum. 

The UA is one of these semi-analytical approaches, that involves using a single approx¬ 
imating function to find the solution of a given differential equation. Although this method 
is quite well known in mathematics as part of the broad area called asymptotic analysis 
[17, 18], it appears to have only found its way into solving the Mukhanov-Sasaki equation 
in 2002 through the work in Ref. [22]. The authors mainly dealt with demonstrating the 
method in calculating the power spectra and spectral indices for the slow-roll single-field 
inflation models and its variants. Specifically, they laid down the working equations for these 
quantities to leading order in the UA and highlighted the advantages of the method namely: 
(a) systematically improvable, (b) has definite error bounds, and (c) does not rely on slow- 
roll parameters being constant. It was later extended and given a more detailed account in 
the follow-up paper where the working equations for the power spectra and spectral indices 
were derived up to next-to-leading order in the UA [23]. The natural progression^ after this 
was the calculation of the power spectra in the more general fc-inflation where unlike that 
of the canonical single-field inflation, the speed of sound is not constant. In Ref. [15] the 
power spectra were calculated up to first order in terms of the so-called Hubble and sound 
flow parameters (see Subsec. 2.2 and 3.2). In a follow-up paper, the power spectra were 
computed up to second order in terms of the same parameters [13]. Note that these last two 
calculations were done to the lowest order of precision in the UA, that is, up to leading order, 
but up to second order in the Hubble and sound flow parameters. Up to second order in the 
Hubble and sound flow parameters, the computations for the higher-order corrections in the 
UA for the power spectra were carried out in Ref. [16], superseding the results in Ref. [13]. 

The problem that we consider in this work does not involve higher-order corrections in 
the UA that were considered in Ref. [16]. Restricting to leading order in the UA, our main 
focus is on the issue of calculating physical quantities to any order in terms of the Hubble and 
sound flow parameters. Consider for instance, the power spectrum for scalar perturbations. 
One has Pg^"^ = Ps{l + ■ ■ ■), where Pg^'^ is the exact power spectrum and the dots represent 

Along the way, other areas of investigation were pursued by several authors. These include studies involv¬ 
ing multiple turning points [24], gravitational quantum effects modifying the dispersion relation corresponding 
to the Mukhanov-Sasaki equation [14], chaotic inflation [25], etc. 
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higher-order corrections in the UA. The term Pg is the base quantity we are interested in, 
that can be computed in terms of the Hubble and sound flow parameters. Since there is no 
practical limitation as to how far the “dots” can be computed, this base quantity Pg limits 

(e) 

up to what order in terms of the Hubble and sound flow parameters Pg ^ can be calculated. 
Furthermore, any defect in the value of Pg is correspondingly reflected in the approximation 
of the exact power spectrum. Now, based on Ref. [13], one has 




Pg = — lim 

2a^eiu. 




( 1 . 2 ) 


where 



with ei = e, so-renamed to fit the definitions of Hubble and sound flow functions (see Subsec. 

2.2) and r/* is the value of the conformal time at the so-called turning point (see Subsec. 

2.3) . The natural way to compute Pg as was done in Refs. [13, 15], is to expand the set of 
quantities c^, a^, ei, and Vg in the coefficient of the exponential function in (1.2) and the 
integrand in (1.3) about the turning point. Consider for instance, ei: 

7 ? 77 1 77 

ei = ei* - ei*e 2 * In — - e?^e 2 * In — - In^ — H-. (1.4) 

r]* V* ^ V* 

Here, denotes the Hubble flow function (e^) evaluated at the turning point. Similar 
expressions can be derived for other quantities involved in the expression for Pg above (see 
Sec. 3). The characteristic form is a chain of Hubble (e^) and sound {6n) How functions 
evaluated at the turning point—the Hubble and sound flow parameters (e„*,<5,i*)—and an 
increasing power of the logarithm of the conformal time. Note that the expansion above for 
ei and other quantities such as o, Vg, etc, can be carried out indefinitely with respect to 
(en*)<^n*) (see Subsec. 3.2). One then performs the integration in (1.3) and long algebraic 
manipulations follow in an attempt to find Pg to a desired order with respect to {en*-,^n*)- 
As already mentioned, in Ref. [15] Pg was calculated up to first order while in Ref. 
[13], it was calculated up to second order with respect to (en*,(5n*)-^ We point out that 
these calculations rely heavily on the striking balance of the logarithmic terms, that is, those 
involving ln(r//ry*), in the exponential term and in its literal coefficients in the expression for 
Pg in (1.2). Indeed, the fact that these logarithmic terms completely cancel out for the first- 
and second-order calculations allows one to take the limit as rj ^ 0~ in (1.2) and finally 
come up with a finite result for Pg. In Ref. [13], such an exact cancellation is considered 
as a “consistency check” of the method that they used because the power spectra involve 
constant modes. This seems to give us no reason to suspect that the calculation of the power 
spectra and other physical quantities such as spectral index and tensor-to-scalar ratio can be 
carried out indefinitely to any order with respect to (e„*, 6n*)- 

Our aim in this paper is to investigate how far one can go with such a calculation 
laid down in Ref. [13] before possibly hitting a divergence, thereby, check its consistency 
when applied to the calculation of the higher order parts of several physical quantities with 

^In Ref. [16], the power spectra were also calculated up to second order with respect to but all 

other quantities including spectral indices, running, and tensor-to-scalar ratio were computed at least to third 
order in 
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respect to (e^*, 6n*), and characterise the nature of the (possible) divergence(s). In particular, 
we consider the calculation of the power spectra, spectral indices, running of the spectral 
index for scalar and tensor perturbations, and tensor-to-scalar ratio, using the UA and the 
expansion scheme employed in Ref. [13]^. We wish to find out up to what order in terms of 
one can calculate these quantities without encountering a logarithmic divergence 
in the end result. We note at this point that although current experiments may not call 
for predictions beyond the second order with respect to (e„*,5n*), our aim here is more 
geared towards mathematical consistency. The hope is to point out limitations in the current 
computational scheme to give way for possible corrections and self-consistent predictions. 

This paper is organised as follows. In the next section, we discuss the set up of this 
work. Specifically, we cover the type of inflation we are interested in namely, fc-inflation, and 
then discuss the Hubble and sound flow functions and the UA. At the end of this section, we 
state the expression for the power spectra for the scalar and tensor perturbations to lowest 
order of precision in the UA. In Sec. 3, we derive a series expansion of the conformal time 
in terms of the Hubble and sound flow functions. The existence of such a series as we shall 
see, plays a significant role in calculating the series expansion of the Hubble and sound flow 
functions with respect to the conformal time and their values at the turning point. With 
the working equations for the power spectra and the expansion for the Hubble and sound 
flow functions in hand, in Sec. 4 we push the limit of the currently available results for the 
power spectra and identify at what order the logarithmic divergences start to appear. Then 
in Sec. 5 we consider other physical quantities such as the spectral index, running, and the 
tensor-to-scalar ratio. Having laid down all the necessary expressions for physical quantities 
such as the power spectra, tensor-to-scalar ratio, etc., we make detailed remarks in Sec. 6 on 
the origin of logarithmic divergences and their ramifications in light of existing calculations 
in the literature. In the last section, we state our conclusion and future prospects in relation 
to the possibility of curing the logarithmic divergences that we encounter in this work. 


2 Set up 
2.1 /c-inflation 

We consider a scalar field 4> called the inflaton, that is minimally coupled to the background 
geometry described by Einstein gravity. The most general local action involving an arbitrary 
function P{X,(j)) where X = can be written as 

S=\J V^d‘^x[M^,R + 2Picj,,X)]. ( 2 . 1 ) 

Here, R is the Ricci scalar, g is the determinant of the metric g^u, and Mpi is the re¬ 
duced Planck mass that we hereafter, take (the square) as unity for brevity. Such an action 
corresponds to the so-called k-inflation. [19, 20] The scope of this work involves slow-roll 
fe-inflation. When P = X — V{4>), where U(</>) is the inflaton potential, the A:-inflationary 
action above reduces to that of the canonical single field inflation. 

The background geometry corresponding to the action above is described by the con¬ 
formally flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric given by 

ds^ = dt^ — a^{t)6ijdx^dxP 

®An equivalent expansion scheme was also used in the calculations in Ref. [16]. 
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With this metric and the equation for the action given above, the Friedmann and continuity 
equations take the form 




SttG 


E, 


E = -3H (E + P), 


(2.3) 


respectively, where E is the energy density, H = d/a is the Hubble constant, and an overdot 
denotes differentiation with respect to the coordinate time t. 

To study the perturbations about our otherwise, conformally flat background, we must 
choose appropriate slicing and threading of spacetime. In connection to this, we find it con¬ 
venient to use the Arnowitt-Deser-Misner (ADM) formalism [26] and decompose the (per¬ 
turbed) metric as 

ds^ = N^dt'^ - hijidx^ + N^dt){dx^ + Wdt), (2.4) 

where N and are the lapse and shift functions respectively. Furthermore, we choose 
the uniform inflaton gauge (also called the comoving gauge) where the inflaton fluctuation 
vanishes, i.e., 5(f = 0, and the spatial components of the metric take the form 

9ij = = 0). (2.5) 

Here, 7^ and are the gauge-invariant scalar and tensor perturbations respectively. The 
tensor perturbations satisfy the auxiliary conditions that they be both transverse and trace¬ 
less. 

It is now possible to write down the equations of motion for the perturbations 7^ and 
For the scalar perturbation, one decomposes the action as S = -h + • • • j where 

5 '(’T-) means the nth-order action with respect to 7Z (see e.g.. Ref. [27]). The background 
corresponds to the zeroth-order part while the interaction starts at n = 3. It is from the 
free part that we derive our equations of motion. In Fourier space, we find the so-called 
Mukhanov-Sasaki equation given as 

'^k + = 0 , ( 2 - 6 ) 


where Vk is related to TZk by Vk = zTZk, the quantity k is the wavenumber, = 2a'^elc^ with 
e = —H/H"^ being the hrst slow-roll parameter, and the symbol prime indicates differentiation 
with respect to the conformal time rj with drj = dt/a. The speed of sound is defined as 


2 Ex 


Ex 


E 


,x 


Px + 2AP 




(2.7) 


where the subscript X denotes differentiation with respect to X. The second equality follows 
from the relation E = 2XPx — P which is due to 00-component of the energy-momentum 
tensor (T^o) being identified as the energy density. For stable solutions wherein > 0, 
corresponding to a bounded total Hamiltonian involving both perturbations and background, 
Px > 0 and 2APxx + Px > 0 (see Refs. [19, 20, 28] for more details). Note in (2.7) that 
when P = X — V we have Pxx = 0 giving = 1, and we recover the usual Mukhanov- 
Sasaki equation for the canonical single-field inflation. For the tensor perturbation, following 
a similar route, we find 


9kE ( k — ) — 0. 


( 2 . 8 ) 
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where /U^/a is the amplitude of gravitational waves. 

Given the two equations of motion above, the (dimensionless) power spectra for scalar 
and tensor perturbations can be written as 


p(e) 

s 


27r2 z 



2fc^ ^ 2 
7r2 a 


(2.9) 


where (e) means exact. It then remains to find the respective solutions (or approximate 
solutions) of the two differential equations above for Vk and //fc to determine Ps ~ Pj^^ and 

(e) 

Pt ~ Pt ^ where Pg and Pt represent approximations based on semi-analytical method such 
as the UA. 


2.2 Hubble and Sound Flow Functions 

When = 1, a closed-form solution exists for the power law inflation [29]. In such a case, 
e = constant and the equations of motion for scalar and tensor perturbations take exactly the 
same form. In general, however, the two equations of motion above do not have a closed-form 
solution.'^ Nevertheless, the structure of the Mukhanov-Sasaki equation makes it amenable to 
semi-analytical calculation such as the UA—the method we employ in this work (see Subsec. 
2.3). 

To proceed with the calculation, one introduces a set of functions or parameters in 
terms of which, physical quantities like power spectra may be expressed. The most commonly 
used is the set of slow-roll parameters defined (albeit, slightly different from each other) in 
Refs. [31-34]. In Ref. [35] a new set of functions called horizon (or Hubble) flow functions 
was introduced. As the authors pointed out, it offers the advantages of being (a) model- 
independent, (b) concise, and (c) easily memorised. All the slow-roll parameters or functions 
defined in Refs. [31-34] can be written in terms of these new functions. Additional flow 
functions called the sound flow functions were introduced in Ref. [15] to include the speed 
of sound in the calculation of the power spectra for /c-inflation. In this work, owing to the 
advantages mentioned above, we make use of these two sets of functions. 

The dimensionless Hubble flow functions {cn) and sound flow functions {dn) are given 
by the recursive definitions [13, 35] 


^n+l 

dn+1 


din €n 

dN 

dlndn 

dA^ 


eo 


H 



( 2 . 10 ) 

( 2 . 11 ) 


where N = f da/a is the number of e-folds, the number n is a non-negative integer, and the 
subscript i indicates some initial value. Note that with the above definitions, ei is simply the 
usual slow-roll parameter e = Furthermore, di serves as a measure of the relative 

rate of change of the speed of sound; ie., di = — (l/P)(cs/cs). 

With the Hubble and sound flow functions in hand, it is now straightforward to rewrite 
the Mukhanov-Sasaki equations for scalar and tensor perturbations. The only relevant parts 

short review of the known analytical solutions can be found in Ref. [30]. The authors of this work also 
presented their own additional analytical solutions. 
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that need modification are the terms z" jz and a!'ja. We identify these as the effective 
potentials Ug and 

z" a!' 

Ug = -, Ut = ( 2 . 12 ) 

z a 

in analogy with the equation for a parametric oscillator. With these identifications and the 
recursive definitions above for and 5n we find 


Ug{v) 


{aHf 


ei + 2 £2 + 4 £2 


^€162 + ^ 6263 + (3 — ei + 62)51 + 5 i + 5 i 52 , 


Utiv) = {aHf {2 - ei) 


(2.13) 


These two equations for the effective potentials are exact. The structure of the expressions 
for the potentials tells us that for the usual slow-roll inflation where f = 1, owing to the fact 
that r] = —l/{aH){l + ei + • • •), the solutions of the Mukhanov-Sasaki equations can be well 
approximated by Hankel functions [31] when the slow-roll parameters are nearly constant. 
Beyond this is a generalization to the case where the slow-roll parameters are not constant, 
inclusion of higher-order corrections in terms of these parameters, and a consideration of 
non-constant speed of sound. The important corrections may be captured by semi-analytical 
approaches utilizing Green’s function [32, 36], WKB approximation [37-39], the UA [13, 16, 
23], etc. 


2.3 Uniform Approximation and Power Spectra 

From (1.1) and the first of (2.12), the Mukhanov-Sasaki equation for scalar perturbations in 
the canonical single-field inflation takes the form 

< + (fc2 - Ug) Vk = 0; Ug = z"/z (2.14) 

When the potential term is negligible in comparison to the wavenumber, the solution of this 
differential equation is simply a sum of two sinusoids. On the other hand, for long-wavelength 
modes, the differential equation approaches a form given by v'^ — {z”jz)vk = 0. It has the 
obvious solution ~ z which corresponds to the frozen modes TZ^ ~ constant involved in 
the power spectrum. More formally, for subhorizon modes wherein |A;/(ai7)| ~ [/cr/] ^ 1, we 
have a sinusoidal solution and for superhorizon modes wherein \k/{aH)\ ~ \k'q\ <C 1, we have 
Vk ^ z corresponding to TZk ~ constant. 

In the case of the more general /c-inflation—the type of inflation we are interested here— 
the effective potential depends on a chain of Hubble and sound flow functions (cf. (2.13)) 
and as pointed out earlier, the resulting differential equation may not have a simple closed- 
form solution. However, the behaviour of the solution at extremal ends of the spectrum of 
krj ~ k/{aH) remains the same; that is, we have two solutions corresponding to subhorizon 
and superhorizon modes and these two are connected somewhere between the two extreme 
values of krj. In principle, one can do a sort of matching typical of a square-well potential 
problems in one-dimensional quantum mechanics. In contrast, such a sort of matching dif¬ 
ferent solutions for different regions is not explicitly existent in the UA. Instead, there is a 
single approximating solution that takes on all the burden of adjustment to fit the boundary 
or initial conditions and the differential equation that it satisfies. 

In the case of the Mukhanov-Sasaki equation, the “differing” behaviour of the solu¬ 
tion at extremal ends in the range of values of krj is taken into account by considering the 
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turning-point. More specifically, one can rewrite the Mukhanov-Sasaki equation for scalar 
perturbation in the canonical single-field inflation as^ 


where 


+ 


1 

Ar] 


-712-9siv) 


= 0 , 


9s{v) = 




(2.15) 


(2.16) 


The turning point corresponds to the value r/ = ??* at which gs{9*) = 0 and the solution 
is in transition from one sort of behaviour to another; say, from oscillatory to decaying. 
Fortunately, owing to the range of values of rj being non-positive, the Mukhanov-Sasaki 
equation has only one turning point. With this feature, it turns out that the solution to 
leading order in the UA can be written in terms of Airy functions. More specifically, one 
takes 


Vk 


lim 
fcrj—>0“ 


01 + 


Bi(/) 


(2.17) 


where ci and C 2 are constants and 


fiv) 


\v-9* 

V -v* 


Q rv ,_ 

- / dry/^ 


(2.18) 


As one can see, in the UA, the problem of solving the Mukhanov-Sasaki equation is converted 
to solving integrals that serve as arguments of the Airy functions in (2.17). The ability to 
find Vk in closed form to some orders with respect to {cn*, Sn*) then relies on our being able 
to calculate those integrals in closed form.® 

In the calculation of power spectra and other physical quantities like spectral tilt 
and tensor-to-scalar ratio, we are only interested in superhorizon modes so we take Vk ~ 
liijikrj^o-if/g)^Bi{f). For large / corresponding to kg —0“, we have the asymptotic form 

Bi(/) ~ e3-^^/(YTr/4). It follows that 


Vk ~ lim 

kr]—¥0~ 



(2.19) 


The generalisation to /c-inflation is quite straightforward. Effectively, one simply re¬ 
places the wavenumber k hy kcg so that the function g is modified as 


9s{g) 


^-k 




( 2 . 20 ) 


The function / and Vk given by (2.18) and (2.19) respectively, take the same form except 
that they involve the more general expression for gg given above. Assuming that the initial 

®Analogous expressions hold for the tensor perturbations. 

®The error corrections in UA also take an integral form; see for instance, Refs. [16, 23]. 
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state of the perturbations is the Bunch-Davies vacuum as in Ref. [13], the power spectrum 
given in the first of (2.9) to leading order in the UA can be rewritten as 





( 2 . 21 ) 


where 



( 2 . 22 ) 


For the tensor perturbations, the power spectrum given by the second of (2.9) can be ap¬ 
proximated as 



(2.23) 


where 



(2.24) 


Note that the limit ry —)• 0 is implicit in the expressions for Pg and Pt above. These two 
expressions for the power spectra are our working equations to leading order in the UA. 


3 Series Expansion of the Conformal time and the Hubble and Sound Flow 
Functions 

3.1 Conformal Time 

With the working equation for Pg and Pt given in the previous section, the possibility of 
calculating power spectra to any order in the Hubble and sound flow parameters mainly 
relies on the convenience provided by the recursive definitions of and 5n given by (2.10) 
and (2.11) and the ability to express the conformal time in terms of On up to any order with 
respect to On- It is often the case that recursive integration is used to expressed rj up to second 
or third order in with the possibility of extending it indefinitely. In this subsection, we 
give a short and simple derivation of a series expansion for the conformal time by indefinitely 
extending the recursive integration by parts usually used in the literature. 

We start with the definition dry = dt/a. With ?y in the range (—oo,0], this definition 
implies that 



(3.1) 


The quantity dt/a = da/{a?‘H) appears again in the integrand in the right hand side so 
following the first one, we perform another round of integration by parts. 



(3.2) 
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Continuing, we find 


where 


OO 


r] = Y.Tn 

n=0 


1 

a 






(3.3) 


(3.4) 


Equation (3.3) is our sought-for series expansion of the conformal time and our main result 
in this subsection. It is essentially a translation from the relatively more incommodious 
language of integrals to a convenient language of differentiation. To verify its validity, one 
can simply differentiate rj with respect to t and show that dr]/dt = Ija in accord with the 
definition dr/ = dt/a. 

The infinite sum above for the conformal time is a Neumann series [40] of the operator 
H~^d/dt acting on 


v = — 

a 




H 


-1 


a\ dt 


-1 


H-\ 


As such, had we started with 


— [l — H ^^^ar] = H ^ 4^ dr/= —, {a ^ 0, H ^ / 0) 


(3.5) 


(3.6) 


we could have ended with the same expansion for r/ given by (3.3) without going through 
recursive integration. 

Using the series expansion for r/ given above and the recursive definitions for the Hubble 
flow functions we find 


To = 1 
fi = ei 

f2 = el + eie2 

^3 = ^1 + 3 e^e 2 + eie| + 616263 

f 4 = ef + 66^62 + 76 ^ 6 ^ + 6162 + 46^6263 + 3616263 + 6 i 626| + 61626364 

! (3.7) 

where fn = —Tn{aH) is 0{e"'). Because of the structure of the recursive definition of 6^, the 
number of terms involved in fn grows as 2”“^. With the above expressions for fn, up to third 
order, the conformal time can be written as 

V = + 6i (l + 6i + 62 + 6^ + 36162 + 62 + 6263) + . (3-8) 

We will use this expression later in an attempt to calculate the power spectrum beyond the 
second order in {en*,6n*)- 
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3.2 Expansion of the Hnbble and Sonnd Flow Functions 

Reference [13] laid down an ingenious expansion scheme to express the quantities involved 
in the power spectra given by (2.21) and (2.23) namely, a, v, c^, e^, and 5n, in terms of r/ 
and (e„*, 5^*)- Since this expansion scheme plays a significant role in the calculation of the 
power spectrum among others, we present here a thorough review based on the discussion 
in the mentioned reference. Along the way, we also demonstrate how such a scheme can be 
used to expressed all the mentioned five quantities above to any order'. 

The two basic ingredients that we need are (a) the series expansion of the conformal 
time derived in the previous subsection and (b) the recursive definition of the Hubble and 
sound flow functions. Proceeding with ei, we start with the relation ii = Hei €2 taken from 
the definition (2.10). Because dr] = dt/a and r] is as given by expansion (3.8) we find 

^ = -ieie2 + 0(e3) (3.9) 

dr] r] 

Integrating from the turning point at ry = ?]* to some arbitrary value of ry gives 

Cl = Cl* ~ ci*e2*T + 0{e\), (3.10) 


where L = ln(7y/ry*). 

To go beyond that given by the last equation above for ei, we perform an expansion 
about the number of e-folds evaluated at the turning point. More specifically, we have 


Cl 


— Cl* + 


dci 

dN 


(AiV) + 


1 d^ei 


(AiV )2 + ... , 


(3.11) 


where AA^ = N — N^. Because d^ei/dN"^ is (n + l)th order, the expansion above for ei 
tells us that we only need an expression for /S.N up to (n — 2)th order in order to find an 
expression for ei to nth order. Needless to say, we only need to find /S.N to first order to 
derive an expansion for ei one order higher than that of (3.10). 

The preliminary expansion for ei given by (3.10) together with the equation for ry given 
by (3.8) is more than enough to express AN to hrst order. In fact, we can go up to second 
order. To see this, we note from (3.8) that 

- = --[l + €l+ef + £i« + 0(e=>)]. (3.12) 

Because dN = da/a, integration yields 


AN — — (l + €1* + Ci* + ci*e2*)T + ^ ci*e2*T^ + ©(e*), (3.13) 


which is good up to second order as we have just mentioned. When we substitute this in 
(3.11) we find ei up to fourth order. 

Cl = Cl* 11 — e2*T — [ei*e2*T — ^ (e^* + C2*C3*) L^] 

~ [(ci*C2* + ci*e2*)T — (ei*e2*C3* + | ei*e|^)L^ 

^From hereon, the word order pertains to order with respect to the Hubble and sonnd flow functions 
(£n,Sn) or parameters {en.»,Sn*) unless otherwise specified, and not with respect to the order of precision of 
the UA. As already mentioned, for the latter, we work only to leading order of precision. 
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+ ^(^ 2 * + + € 2 *^ 3 * “I" +C>(£*)} (3-14) 

Owing to ei being in the denominator of the right hand side of (2.21), this fourth-order 
expression is just sufficient to calculate the power spectrum Pg up to third order. 

Following a similar approach, we can also find 62 up to fourth order, 

62 = 62*{1 ~ ^3*L — [ei*e3*L — -|- e3*e4*)L^] 

~ [(61*63* -|- 6i*e2*63*)L — (^ 61*62*63* -|- 6i*e|* -|- ei*e3*e4*)L^ 

4" g (63* + 3e|* 64*-|-€3*64*-|-63*64*65*) L^] -|-0(6*)}, (3.15) 

then 63 up to fourth order, and so on. Once we have these, we can update our expression 
for AN up to fourth order, thanks to the infinite series for rj. Then our expressions for 6 „’s 
can be updated to higher order. The iterative approach continues and we can express all the 
Hubble flow functions to any desired order in terms of the Hubble flow parameters and the 
conformal time. 

We employ the same method for the expansion of the sound flow functions with respect 
to the Hubble and sound flow parameters and the conformal time. One starts with the 
recursive definition given by (2.11). For <5i this leads to 

5i = 5u - 5 iJ 2*L + 0{6l) (3.16) 


Then we use the series 




(^ 1 * -|- 


d<5i 

w 


{AN) + 


1 d 2 (ji 


(AiV )2 + ... , 


(3.17) 


The same goes for the calculation of other sound flow functions. There is no limit as to the 
order by which we can express 6n in terms of (5„* and 6 ^* since AN can be determined to any 
desired order. For Ji and <52 involved in the effective potential given by the first of (2.13), we 
have up to third order. 


51 = (5i*{l — S2*L — ^ [261*52*71 — (^2* -|- 52*53*)T^] -|- 0((6*, 5*)^)} 

52 = 52*{1 — 53*L — ^ [261*53*71 — (5|* -|- 53*54*)L^] -|- 0((6*, 5*)^)} (3.18) 

Of the five quantities (a, 1 /, c^, 6„ and 5^) we mentioned near the beginning of this 

subsection, we are done with two namely, 6^ and 6n- For the remaining three, we start with 

the scale factor. Since AN = ln(a/a*) then a = a*e^'^. Because the lowest order in the 
equation for AN is —L = — ln(7//r/*) then 

70 pAA^+L 

= + . + ■■■), (3,19) 

rj H^rj ' 

where in the second equality, we have used the equation for the conformal time. Squaring 
this yields 

2[6i* — 6i*71] -|- [36i* -|- 261*62* — (06^* -|- 26l*62*)71 

+ (26 i * -|- 61*62*) 71 ^] -|- [ 46 ^* -|- 86^*62* + 261*6^* -|- 261*62*63* 
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~ + 14e^^e2* + 2 ei*e|^ + 2 ei*e 2 *e 3 *)L 

+ (Sef* + 9e^^e2* + ei^e^* + 


— 5 ( 4 ef^ + 6ef^e2* + €1*62^, + €1*62*63*)L^] + 0((6*, ( 5 *)^) |. 
The next is the (square of the) speed of sound, c^. We have 


2 2 , 


diV 


(AA^) + 


1 d^c? 


2 !dA ^2 


{/\Nf + ■■■ . 


By virtue of the recursive definition for 6n given by (2.11) we find 

d^c? 


d(5o c- c- , dc2 ^ 2 ^ 
— = 5o5i then _ = -2 c,5i, 


+ 2(5i(52). 

It follows that 

C 5 = + 2(5i*T + [2ei*5i*L + (2(5i* — 5i*(52*).b^] 

+ 2(ei^5i* + 6i*e2*(5i*)T — (6i*e2*(Ji* + 2ei*(5i*(l2* ~ 4ei*(5i^)ii^ 


(3.20) 


(3.21) 


(3.22) 


+ g (ddf* — 6(5 i^(52* + <5i*<^ 2* 4“ + Cl((f*) <5*)^) j"- (3.23) 

Finally, for the index function Ug one can perform an expansion analogous to that for (see 
(3.21)). Alternatively, noting that = j + rfUg-, where the product rj^Us by virtue of (3.8) 
and the first of (2.13) is expressible solely in terms of the Hubble and sound flow functions 
without the extra factor (aH), we can simply use the expansions for these functions we have 
derived to write down the needed expression for Both of these two methods give the same 
result: 


= I + [ei* + ^ 62* + ( 5 l*] + [61^ + 6i*(Ii* + ^ €1*62* + ^ 62*63* + ^ ( 5 i* 52 * 

~ (61*62* + ^ 62*63* + (5l*(l2*) a] + [6^^ + ^ 61^62* + ^ 6i*6|^ 

+ ^ 61*62*63* — ^ 62^,63* + 6^*(5 i* + ^ €i*€2*(5i* — g 62*63*(5 i* 

+ I 6 i* 5 i*( 52* — ^ 62*5 i*( 52* ~ I ~ (36i^62* + ^ 61*62* 

+ 5 61*62*63* + ^ 62*€|^ + ^ 62*63*64* + 6i*62*(5i* + 2€l*(5l*(52* 

lrr2 lrrr\r /I 2 1 1 2 

+ 3 0 i* 02 * + 3 01 * 02 * 03 *) L + ( 2 61 * 62 * + 2 C1*£2*63* + 4 62 * 63 ^ 

+ 5 62*63*64* + ^ <5i*(52* + 5 <5i* 52*<^3*).^^^] + O((€*, d*)"^) (3.24) 

Similarly, for the tensor perturbation index function we have 

z^t = I + 61* + [6^^ + I 61*62* — 6 i*62*T] + [64^ + ^ 64^62* + I 61*62* 

+ I 61*62*63* — (361^62* + I 61*62* + I 61*62*63*) L 

+ 61 * 62 * + ^ 61 * 62 * 63 *)+ 0((6*, (5*)^) (3.25) 

With these last two equations settled, we are now in a position to calculate the power spectra 
for scalar and tensor perturbations beyond the second order. 


- 13 - 






4 Logarithmic Divergences in the Power Spectrum 

We first consider the power spectrum for scalar perturbations given by (2.21) 


Ps = - 






Svr^ a^eiK 


(4.1) 


Before we delve into the possible logarithmic divergences that may arise from this expression, 
noting that we have to take the limit —?• 0“, we address the possible problem coming from 
the factor rj in the right hand side. In fact, by virtue of the expansion for given by (3.20) 
we have a total factor of such that the zeroth-order power spectrum is 


PP = 


127r2 


£ll^^3(g2vI/,)(0) 

^1* 


(4.2) 


A factor of l/rj^ should then fall out of the exponential term to avoid a trivial result. That 
this is so the case can be easily seen from the following calculation. 


2T, = 2 






+ ••• 


2^5 = -3 + 31n2 - 3L H-. 


(4.3) 


Exponentiation then yields 

= + ...), (4.4) 

which cancels as advertised. 

We now go to the calculation of Pg. With all the literal coefficients of available 
from the previous section, all that is left is the evaluation of the exponential term. In order 
to do so, one expands the integrand about (c^^, Ug*)- 


9s 



+ 


C-s* 


dvg 


(Aug) + 


del 


(Ac2) + ... . 


(4.5) 


Here, Ac^ = and Az/^ = Ug — Vg^. Then can be rewritten as 

dutt 


where 


4' ,8 = -i^g 


f^du 2 \i ^s* duu 

J —(1-u2)2+^/ ^ 

r du \ f 

/ / 

Ji u(l — vd) 2 ° Ji 


(1 — v?)^ 

du 


3 (A«)^ 


7 /( 1 — 7X^)2 o Jl (l — u ‘^)2 

duu , . . ,. . duu^ 


1 r 3 (a»)(a..)+^ r 

2ji 16 A 


(As)3 + 


(4.6) 


kCoikT" 


u = —- 


and 


(4.7) 
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From (3.23), (As) is at least first order while from (3.24), (Az/^) is at least second order. It 
follows that the remaining terms following the six integrals above are at least fourth order. 
Needless to say, the visible series of integrals in (4.6) is enough to calculate Ps up to third 
order. 

We further note that all the terms contained in As and Avg involve powers of logarithms; 
i.e., L” = liP u, where n > 1. As such, all the integrals above take the form 

duu^ (^1 — u^y liP u, (4.8) 



where p, q, and r are all integers. This form allows us to determine the integrability of all the 
terms involved in the Taylor series expansion given by (4.5) and see where the logarithmically 
divergent terms may come from. The first term in this series expansion, by virtue of (4.3), 
leads to a logarithmically divergent integral. For the remaining terms, we have the following 
classification : those involving (a) (Ac^)” only, (b) {Avg)'^ only, and (c) {Aug)"'{Ac^)^, where 
n, p > 1. 

The integrals involving only (Ac^)"" or (As)”' take the form 



^2n-l ^ 

(l-u2)"-5’ 


m> n 


(4.9) 


The integrand is continuous on (0,1). Furthermore, the limits as 'U —>■ O"*" and u —)• 1“ 
both exist and vanishing. It follows that all integrals involving (As)” only are finite as 
w —>• O"*". Similar analysis indicates that all integrals involving the product (As)”(Ai/s)^, 
where n,p> 1, are also finite. The integrals involving only (Ai^^)” take the form 



Cl In”^ u 

u[i — v^y ^ 


+ 


C2 In"^ u 

u{i — v?y ^ 


+ ••• 


m > n 


(4.10) 


where Cj’s are constants. The number of terms inside the pair of square brackets depends 
on whether n is odd or even. If it is odd, we have (n + l)/2 terms while if it is even, we 
have (n + 2)/2. The integrand is continuous on (0,1) and the limit as tt —1“ is finite 
and vanishing. The limit however, of every term inside the pair of square brackets above as 
tt ^ 0+ diverges. This leads to 


r du u 

1 u(i — v?y 2 


ln’”+^ w 
m + 1 


+ const. 


(4.11) 


in the limit as tc —0"^. Hence, (a) all divergences involved in are logarithmic in nature®, 
and (b) only those integrals involving (Az/^)” and not (As)” give rise to such logarithmic 
divergences. 

Had we been only interested in logarithmic divergences in Pg, we could have simply 
expressed c^, a?, Vg, and ei in exponential form and combine it with so as to simply 
find whether we could have incomplete cancellation of terms involving Inzc without having 
to compute the finite integrals in (4.6). Nonetheless, for completeness, we compute the full 

®Intermediate calculational results may involve polylogarithms of the form Lisfconst/u;") which in the 
limit as w —>■ 0^, result to logarithms. 
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power spectrum at least up to third order. In order to do so, we perform the necessary 
integrations in (4.6). Here are our results after a lengthy calculation. 


Aq = 

pw 

1 

Bi = 

nw 

1 

B 2 = 

pw 

1 

B 3 = 

pw 

1 

Cl = 

pw 

1 

C 2 = 

pw 

1 

D 2 = 

pw 

1 

Ds = 

pw 

1 

E 2 = 

pw 

1 

F 3 = 

pw 

1 


u 

u \nu 

(1 — u2)2 


= 1 — In 2 + In re 


= 1 - In 2 


TT 


-T = -2 - In^ 2 + 2 In 2 + 

(l-a2)i 12 

=-Ic(3)-61n2 + 6+^ln2 + 31n^2-lnl2-^ 

{l-y2)-2 2 4 4 


Inu 


tt (1 — u2)2 


= — -- In^ 2 + - In^ re 

24 2 2 


In^ u 


u{l — v?)'^ 
V? \r? u 


1 , , vr^ 1 1 q 1 1 q 

1 =--C(3) + -ln2--ln3 2+-ln3u; 


TT 


(l-n2)2 

V? In^ u 
(1 — u2)2 


, =2 + 21n^2-21n2- 

2^i 6 


2 2 
= 3C(3) - 6 - 3 In^ 2 + 6 In 2 - — In 2 + 2 In^ 2 + ^ 


TT 


, = +ln2 2 

/2if 12 


^ = -4C(3) +6 + 21n22 - -ln3 2+ — ln2 - 61n2 - — (4.12) 

(l-n2)2 3 3 6 

The labels A, B, ■ ■ ■ ,F correspond respectively, to the visible six integrals in the right hand 
side of (4.6). Moreover, the subscripts of these labels indicate the power of the log function 
in the integrand. 

Using the equations for definite integrals above in (4.6) and then substituting the result¬ 
ing equation together with the expansions for ei, a^, c^, and Us given by (3.14), (3.20), (3.23), 
and (3.24) respectively, in the working equation for the power spectrum given by (4.1), and 
noting that kcg* = —i's*/v* 1^6 turning point, we find 


where cq = 18e ^ 
= 1 


“ 87r2 


0.896 and 


H^co 




.[p]0)^p(l)^p(2)^p(3)]^ 


(4.13) 


= - (I - 2 In 2) ei* - (g — In 2) € 2 * + (| - In 2) di* 
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^ ln2 + 21n^ “ (t “ “ 5 2 )ei*e 2 * 

- (^ + 5 i \v? 2 )el^ - (i - - i ln2 + i In2 2)e2*e3* 

- (f - f \n 2 + 2\n^2)eu5u - (f - f In2 + In^ 2 ) 62 *^ 1 * 

+ (f| - I ln2 + i In^ 2)4 + (f - ^ vr^ - I ln2 + i In^ 2 ) 5 iJ 2 * 

= (I - I In2 - I In^ 2 + I In^ 2 )ef^ + (I + ^ vr^ + i TT^ In2 - 6In2 + 3In^ 2) e\^e 2 * 

- (f -|^^-4(3) + iln2-4n22-lln32)6i*ei 
+ (^ - 11 1^2- I ln22 + i ln3 2)ei 

- (^ - vr^ - 4(3) - 21n2 + I In^ 2 + I In^ 2 )ei*e 2 *e 3 * 

- vr^ - ^ vr^ ln2 + 4^2 + 4n^ 2 - 4n^ 2)e2*e3* 

{72 + i C(3) — ^ vr^ ln2 — I In^ ^ ^ 2 )e 2 *e|^ 

+ (4^^ + i ^(3) - ^7r^ln2 - i In^ 2 + 4n^ 2 ) 62 * 63 * 64 * 

- (I-21n22-I In2 + 21n3 2)6?,(5i* 

- (I + ^ TT^ - ^ ln2 + TT^ ln2 + 4n^ 2 + In^ 2 ) 6 i* 62 *( 5 i* 

“ (I “ I l’^2 - 4n^2 + i ln^2)6i^5i* 

+ (In^ 2 - In^ 2 + i - i ln2)6i*4 + (| “ 4^2 - 2 + 2 ) 62*4 

+ TT^ + 4^2 - ^ TT^ ln2 - 4n^ 2 + 4n^ 2 ) 62 * 63 *( 5 i* 

“ (^ “ 11 2 — 4^^ 2 + 4n^ 2 ) (5f^ 

— TT^ — ^ ln2 + ^ TT^ ln2 + 41n^ 2 — In^ 2 )eu 6 u 62 * 

- vr^ — 4^2 + ^ vr^ ln2 + 21n^ 2 — 4n^ 2 ) 62 * 6 ^ 62 * 

+ TT^ - 4^2 + ^ TT^ ln2 + 21n^ 2 - 4n^ 2)4<^2* 

+ (3 - - 4(3) - 31n2 + ^7r2ln2 + 4n^2 - I ln^2)(5i*4 

+ (3 - ^ TT^ - 2 C(3) - 31n2 + ^ TT^ ln2 + g In^ 2 - g In^ 2)5i*(52*53* 

+ 42 ^ 1 *^ 2 * + 261 * 62 * 63 * + 62 * 6 ^^ + 62 * 63 * 64 * + 25i*4 + 25i*( 52*(53*) Inw. (4.14) 

The result above for the power spectrum generalises the second-order approximation for Pg 
given in Ref. [13]. The third-order part of (4.13) is something new. What is surprising about 
this part (see the expression for Pg '), is the unexpected set of terms involving Inic which is 
divergent in the limit as rc ^ 0'*'. Unlike that of Pj^\ Pi^\ and pj^\ we have an incomplete 
cancellation of (divergent) log terms. Furthermore, considering the case of a constant speed 
of sound does not help eliminate these terms. Since terms involving Inrc are carried to the 
calculation of the higher-order power spectrum, we expect Pg ^ for n > 4 to also contain In w 


17 


(raised to some positive integer). 

The calculation for the power spectrum for the tensor perturbations is less cumbersome 
than that of the scalar perturbations. The speed of sound is unity and out of the 10 integrals 
used for Pg, we only need three corresponding to the first and third terms in the right hand 
side of (4.6). Following a similar approach, now with (2.23) as the working equation for Pt 
instead of (4.1) for Pg, we find 


Pt = 



p(0) 





(4.15) 


where 


p{0) 


= 1 


= (-| + 21n2)ep 

= (f - f In2 + 21n2 2)e?, - (f - | In2 + In^ 2)epe2* 

Pf = (|-|ln2-|ln22 + |ln32)eL 

+ (^ TT^ — ^ ln2 + 4 TT^ ln2 + 51n^ 2 — 21n^ 2 )e^^e 2 * 

— (t “ i ^ “I ln2 + ^ vr^ ln2 + I In^ 2 — | In^ 2 )ei*e 2 * 

— (t “ i ^ “ i ln2 + ^ vr^ ln2 + I In^ 2 - 4 In^ 2 )ei*e 2 *e 3 * 

+ + ei*e2*e3*) Inrc (4.16) 


Just like that of the power spectrum for scalar perturbations we have an incomplete cancel¬ 
lation of (divergent) log terms. 

Let us try to look further into the details of the emergence of these logarithmically 
divergent terms in the power spectrum for the scalar perturbations, if only to convince us of 
such incomplete cancellation. (A similar analysis holds for tensor perturbations). In (4.14) 
there are six terms involving Inrc. 

4) (|ci*e2* P + ^£2*63* + ^£ 2 *£ 3 *£ 4 * + ’^^1*^2* + |<^i*<52* J3*) ln(4-17) 


Terms in the working equation for Pg contributing to this set of log terms may be seen by 
rewriting (4.1) in exponential form as alluded to earlier. To do this, use the dehnitions of the 
Hubble and sound flow functions for ei and c^. For a? note that da/a = — (d??/r/)(l + a) where 
the conformal time is written as r] = —(1 + a)/{aH). For Vg observe that l\Vg = Vg — 
is at least second order so (Az^^)^ can be truncated in the exponentiation. With these 
considerations, we can rewrite ei, c^, a?, and Vg in exponential form as 


ei = ei* exp 


cl = cl^ exp 


r dr 

- — £2(1 + a ) 

Jri, r 

n dr 

- 2 / 
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exp 


-2 


2 ^ (1 + a*)"^ 

1 1 _A. 

— = —e 
1/ 



(4.18) 


where the last equation only holds up to third order. Substituting into the working equation 
for Pg given by (4.1) we find 


where 


Ps = N{r]) exp 



a{2 + €2 + 26i) + 62 + 25i 


Az/ 

-h 2'I'5 (, 


Af = 


(1 + ct*)^ 


(4.19) 


(4.20) 


We can now focus on the three-term argument of the exponential function in (4.19). 
Specifically, noting that all third-order terms involving In^ w and In^ w are canceled along 
the way, we pay attention to those third-order terms involving only the first power of Intc. 
Considering 24^5, a third-order term involving the first power of In re can only come from the 
first integral in the right hand side of (4.6) since all integrals involving (As)"' are convergent 
and the divergent integrals involving (Az/^)" start at In^ w (see the fifth and sixth of (4.12)). 
It follows that 2'ks only contributes the following third-order logarithmically divergent terms 
involving the first power of Inzu: 


9U/ ^ (‘->P a a a _L 34 ^ ^2 , 77 ^2 ^ 

3 — I ^6]^^ 61*62*63* -g 61*62* -r 61*62* 

~ 5 £ 2 *^ 3 * + 26i*hl* -|- I 6l*hl*(52* + ^ 6 i*62*5i* 

- I 62*(5l*h2* - I 62*63*hi* - | (5i*h2*) Ill'll (4.21) 

The second term in the argument of the exponential function in (4.19) contributes 

^/l4 2 I 2,4 12 ,2 IT 

-D I -g- 6 ]^^ 62 * -h 61 * 62 * -h 3 61 * 62 * 63 * — g 62 * 63 * -h g 61 * 62 * 01 * 

+ I 6l*hi*52* — I 62 * 63 *^ 1 * — I 62*(5i*52* — I 

+ 5 ^ 2 * 03 * + ^ 62 * 63 * 64 * -|- I hi*52*<^3* + I (4.22) 

Already at this point, one sees the last four terms above inside the pair of parentheses match 
the last four terms inside the pair of parentheses in the expression for Pg given by (4.17). 
These four terms are nowhere to be found in the contribution of 2Ts. Furthermore, writing 
out the integrand in (4.19) as 

0(2 -|- 62 -|- 2 ( 5 i) -|- 62 -|- 2 hi = ( 26 i -|- 62 -|- 25 i ) P (86162 -|- 26i -|- 2 ei 5 i ) 

P (2616263 -|- 86162 -|- 7efe2 T 26^ -|- 26i62(5i -|- 263 ^ 1 ), (4.28) 


tells us that with the current expansion scheme, there is no way to produce the (negative of 
the) mentioned four terms out of the integration of the last equation above. For instance, 
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5i*( 52*(53* In w can only possibly be produced from the third term (2(5i) in the right hand side of 
(4.23). But the expansion of 26i only leads to the closest possible term In^ w which 

upon integration yields | 5i*(52*53* In^ w. Upon further considering the remaining three terms 
namely, ^ g e 2 *e 3 *e 4 *) and | we find that the second argument of the exponential 

function in (4.19) namely, —{Ai's)/i's, is the sole contributor of four out of six logarithmically 
divergent terms in Pg . The remaining two are a product of incomplete cancellation from 
the three-term argument of the exponential function in (4.19). 

Hence, with the expansion scheme laid down in Ref. [13] (thoroughly reviewed and 
elaborated in Sec. 3.2) and the method of the UA to leading order of precision, the power 
spectra for scalar and tensor perturbations fail to converge starting at the third order with 
respect to the Hubble and sound flow parameters. Looking back, we see that the expansion 
scheme already involves divergent terms in the limit as ^ 0~. For instance, 

ei = ei* - ei*e 2 * In — H-. (4.24) 

One then delays taking the limit to the end of the calculation of the power spectrum in the 
hope that all the possibly logarithmically divergent terms would cancel along the way. It 
turns out that this technique only works up to second order. 

We wish to emphasise that the failure to converge cannot be attributed to the UA itself. 
The working equations for the power spectra given by (2.21) and (2.23)— the main machinery 
of the UA in calculating Pg and Pt — are mathematically well-defined (finite) expressions. 
This leads us to put the blame on the implementation of the working equation itself. When 
applied to the case of the slow-roll ^;:-inflationary scenario to obtain quantities appropriately 
decomposed with respect to the Hubble and sound flow parameters, this implementation 
requires a suitable expansion scheme for the slow-roll parameter, speed of sound, etc. As 
we have found, this scheme has brought with itself terms that can blow up in the limit as 
r] ^ 0~, that eventually fail to completely cancel in the end result starting at the third order 
for the power spectra. To be clear, the origin of the divergences is the expansion scheme 
itself, not the UA. 

It is worth noting that an “alternative” expansion scheme was considered in Refs. [14, 
41]. It involves series expansion about the turning point but the variable of differentiation is 
the conformal time instead of the number of e-folds in the current scheme. For instance. 


= l^g* + 


di^g 

dry 


{r]-V*) + 


1 d'^i^g 
2! dry2 


iv - + ■■■ 


(4.25) 


However, (although the proof requires a lengthy calculation,) this expansion scheme is equiv¬ 
alent to that covered in Subsec. 3.2; therefore, it is not an “alternative” one. Furthermore, 
using this in its current form is practically cumbersome at least for the inflation model we 
are considering here. In agreement with Ref. [14], the higher-order derivative terms for 
slow-roll inflation can contain lower order terms with respect to the Hubble and sound flow 
parameters. It follows that ignoring higher-order derivative terms amounts to a significant 
loss of precision. References [14] and [16] adopted the current expansion scheme (discussed 
in Subsec. 3.2) to avoid this problem in their calculations of the power spectra up to second 
order in (e^*, 5n*) and up to third-order precision in the UA. 
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5 Tensor-to-scalar Ratio, Spectral Index, and Running 


5.1 Tensor-to-scalar Ratio 

The tensor-to-scalar ratio is defined as r = Pt/Pg- Since the logarithmically divergent term 
In re starts to appear in the third-order part of both (approximations for) Pg and Pf, one may 
expect that the third-order part of r should contain In re as well. That this is not the case 
can be seen by substituting the expressions for Pg and Pt given by (2.21) and (2.23) in the 
definition of r above. 


r = 16%—(5.1) 
cj yt 

Here, the limit as ?? —)• 0“ is again implicitly assumed. The term ei in the denominator in the 
expression for Pg is brought up making the zeroth-order part of r vanishing. As a consequence, 
the tensor-to-scalar ratio can be computed up to third order without encountering logarithmic 
divergence. Here is our result. 

r = 16ei*Cs* -I- -I- -I- , (5.2) 


where 
f(0) = 1 

= (3 - In 2 )e 2 * - (I - ln2)(5i* 

f( 2 ) = (I - In 2 )ei*e 2 * -F (g - 5 ln 2 -h ^ ln^ 2 )e 2 * + [\ - ^'k'^ - \ \n2 + \ In^ 2 ) 62 * 63 * 

~ (t ~ 31n2)6i*(5i* — (I — ^ ln2 -I- In^ 2 ) 62 *(^i* -|- (^ — ^ ln2 -|- | In^ 2) <5^^ 

- (^ - ^ TT^ - I ln 2 -F i In^ 2 )( 5 i*( 52 * 

,-(3) = (25 _ ln 2 ) 6 L 62 * + (i - 2 ^ vr^ - f ln2 + I In^ 2 ) 61 * 6 ^ 

+ (m “ S + S “ I ln^2)6^* + “ I ln2 -Fln^ 2 ) 61 * 62 * 63 * 

+ (^ “ 7 ^ “ 5 2 -F ^ vr^ In 2 -b ^ In^ 2 - i In^ 2 ) 63 * 63 * 

- (^vr^ + iC(3) - ^vr2ln2 - i In^ 2 + i ln3 2 ) 62 * 6 i^ 

- vr^ + 5C(3) - ^ 7r^ln2 - I In^ 2 + ^ In^ 2 ) 62 * 63 * 64 * 

- -51n2)6U*+(f-ii2ln2 + 31n22)6i*<^L 

- (^ - ^ TT^ - 16 In 2 + ^ In^ 2 ) 6 i* 62 *(^i* 

+ (§ - f ln2 + I ln2 2 - i ln3 2 ) 62 * 5 ?* - (^ - ^ In2 + f In^ 2 - 1 In^ 2)6i,5i* 

- ^ TT^ - ^ In 2 + ^ TT^ In 2 + I In^ 2 - i In^ 2 ) 62 * 63 *^ 1 * 

- (^ - i vr^ - f In 2 + 2 In^ 2) 6i*5i*52* 
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- (if “ 11 “ T In 2 + ^ TT^ In 2 + 3 In^ 2 - i In^ 2) € 2 * 6 ^ 82 * 

- ln2+f In22-iln32)4 

+ (^ - I vr^ - In 2 + ^ TT^ In 2 + 5 In^ 2 - i In^ 2 ) 

- (3 - ^ vr^ - 1 C(3) - 31n2 + ^ vr^ ln2 + I In^ 2 - i ln3 2)5u6l, 

- (3 - ^ TT^ - i C(3) - 31n2 + ^ vr^ ln2 + I In^ 2 - i In^ 2)5i*(52*(53* 

“ ^ (c2*C3* + C2*C3*f4* + 2(5 i*(52* + 2(5i* 52*<53*) Inzn. (5-3) 

As can be seen, the presence of terms involving In w in given above makes the fourth-order 
part of r logarithmically divergent. 


5.2 Spectral Index and Rnnning of the Spectral Index 

The spectral index for scalar perturbation defined as 


, _ dlnP, 
dlnA: ’ 

can be computed from the derived expression for the power spectrum in Sec. 4 given by 
(4.13). The idea is that from the relation k = —i's*/i'n*^s*), telling us that ry* is an implicit 
function of k, one can determine dln(— 7 y*)/dln A: (see Refs. [13, 16]). The resulting expression 
can then be used to carry out the calculation of as dictated by the definition above. One 
may think that because the power spectrum under the current scheme involves logarithmic 
divergences starting at the third order, Ug may inherit the same bad feature. However, rj 
and k are independent in the limit as ry —)• 0 “ unless the turning point relating k and ry* 
happens to be at ry = ry* = 0, which is a trivial scenario. Hence, even if we accept that the 
derived power spectrum containing logarithmically divergent parts is valid (which it is not), 
the resulting spectral index will still be finite. 

Let us try to see the finiteness of the spectral indices under the current calculation 
scheme in a more rigorous way. We start with the definition of Ug above with Pg given by 
(4.1) and derive the appropriate working equation. The derivation of the working equation 
was done in Ref. [22] for the case = 1 and we are here going to follow the same path. 
Taking the logarithm of the power spectrum, we find 



dlnPs 
din /c 


3+ 2^:77 / dr yJgg{k,T), 


dk 


which implies upon using the definition of gg given by ( 2 . 20 ), 


(5.5) 


rri Ar 

ng-l = 3-2k‘^ ^ ^ . (5.6) 

4* V9s{k,T) 

This is our sought-for working equation for the spectral index. 

References [22, 23, 41] discussed the square-root singularity at the turning point in the 
integrand above for the case where = 1. For the case in hand, where 7 ^ constant, 
we have the same form of singularity and we wish to find out using the expansion scheme 
used in Ref. [13] if the working equation above for Ug leads to a finite or infinite results 
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with logarithmic divergences just like that of the power spectrum. To proceed, we perform 
a Taylor series expansion in the same way we did in Sec. 4. 


V9s{k,T) 


+ 


d 


Us* ,Cs* 


di^s ^/Fs 


d 


{AUs) + 2 ^ 


(Acf) + 

Us* jCs* 


(5.7) 


The first term in the right hand side corresponds to a convergent integral. The rest of the 
terms can be classified into three types: those involving (a) (Acg)"’ only, (b) (Ai^g)"' only, 
and (c) {Avs)"'{Ac^y, where n, p> 1. 

For terms involving only (Acg)”, with n > 1, the corresponding integral takes the form 



Cl In^ u C 2 liT u 
(1-1x2)’^"^ (1-1x2)’^+^ 


(m > n) 


(5.8) 


where ci and C 2 are constants. The first term inside the pair of square brackets is continuous 
on (0,1). Furthermore, the limits as xx —O"*" and xx —>■ 1“ both exist and vanishing. The 
same goes for the second term except, when m = n, where the limit as xx ^ 1“ diverges. 
This however, does not lead to a divergent integral. That this holds can be seen by using a 
comparison test for improper integrals with the case for n = 1 as the basis. The second term 
can result to a divergent integral if m < n which is clearly not the case here. 

For terms involving only (Azx^)"', with n > 1, the corresponding integral takes the form 



cixxln^'xi C2Xtln™'xx 

(l — xx 2)"'"'~2 (l — xx 2 )"' 2 


(m > n) 


(5.9) 


The number of terms inside the pair of square brackets depends on whether n is odd or even. 
If it is odd the number of terms is (n + l)/2. If it is even the number of terms is (n + 2)/2. 
All the terms inside the pair of square brackets are continuous on (0,1) and the limits as 
XX —7> O"*" and xx —)• 1“ both exist and vanishing except for the first one where the limit as 
XX —)• 1“ diverges. However, as in the case of terms involving (Ac2)”, this does not lead to 
a divergent integral. That this holds can be seen by using a comparison test for improper 
integrals with the case for n = 1 as the basis. The integral can diverge if m < n which is 
clearly not the case here. 

The analysis for cross terms in the Taylor series expansion given by (5.7) involving 
(Axxs)”'(Ac2)^’, with n, p > 1, leads to the same finding that the corresponding integral is 
convergent. This together with the results above for the cases of (Ac^)”" and (Azx^)”, leads 
us to conclude that the special form of the working equation for Ug given by (5.6) enables 
the calculation the spectral index to any order with respect to the Hubble and sound flow 
parameters without encountering logarithmic divergences involving the limit as xc —0"*“ of 
Inxc. 

The examination for the spectral index for tensor perturbation follows the same route. 
From its definition given by 


nt-l = 


din Ft 

dln/s ’ 


(5.10) 


we arrive at the working equation similar to that of (5.6) for the scalar perturbation. 


rax - 1 


= 3-2k‘^ 



dr 

Vgt{k,T) 


(5.11) 
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As in the case of scalar perturbations, the placement of being the denominator of the 
integrand makes the integral above convergent. More formally, the Taylor series expansion of 
the integrand takes the same form as that of the scalar perturbation with = 1 because the 
index functions Ug and vt have the same form. With n* being computable under the current 
scheme up to any order with respect to {en*-, the same holds for rit. 

The fact that and nt can be computed to any order tells us that the running of these 
spectral indices defined as 

_ dn^ _ dn* 

d In A: ’ d In A ’ 

for scalar and tensor perturbations respectively, can as well be calculated to any order in 
(^n*, One method, already elaborated at the beginning of this subsection, is to simply 
use the relation k = —Vs*/{'n*Cs*) for scalar perturbation to find dln(—r/*)/dln A and then 
substitute the result together with the expression for Ug coming from (5.6) in the dehnition 
given above for Og. An equivalent procedure holds for the tensor perturbation. The other 
method involves directly differentiating the expression for Ug given by (5.6) and deriving the 
working equation for ag in the same way we did for Ug. However, the resulting expression for 

3 

this one involves an integrand with gg being part of the denominator different from that of 
1 

Qs in the expression for Ug. This leads to some practical difficulties because of the resulting 
divergent integrals. 


6 Remarks on Logarithmic Divergences 

In what follows we clarify some of the questions or reiterate some of the main points stemming 
from the emergence of logarithmic divergences that we point out in this work. 

The finite lower-order results. Given the exact same method to calculate the first-, 
second-, and so on, order part of a given physical quantity, (in fact, they are all tied together,) 
the emergence of logarithmic divergences say, in the third-order part, questions the validity of 
the lower-order parts even though they converge at second order. As such, the “consistency 
check” in the calculation of the power spectra up to second order carried out in Ref. [13], 
for which all terms that can lead to logarithmic divergence canceled, should be viewed with 
caution; as we have demonstrated, all terms starting at the third order do not pass this 
“consistency check.” It is in this same spirit that the results of calculations carried out in 
Refs. [15] and [16] should be taken with caution. Moreover, it may be worth pointing out 
that even as far as the third order of precision in the UA, the percent difference of some of 
the coefficient of the second-order combinations of Hubble flow parameters, with respect to 
that obtained using other approximation methods, can still be significant; e.g., for ei*e 2 * in 
the expression for the scalar and tensor power spectra, where the symbol * means evaluation 
at the horizon crossing, the percent difference with respect to that obtained through the 
Green’s function method is 90% and 13%, respectively (see Table HI of Ref. [16]). The 
resulting practical expressions for the power spectra and tensor-to-scalar ratio do not rest on 
a solid mathematical foundation and hence, are in need of some further justification. 

We note that this work can only go this far. There is still a substantial possibility that 
the inconsistency in the expansion scheme can only affect some but not all of the finite first- 
and second-order terms for the power spectra and the finite first- to third-order terms in the 
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tensor-to-scalar ratio. A thorough check should be put in place. In order to do this, one may 
need to carry out the calculations well beyond the third order of precision in the UA. The 
succession of the computations however, for the power spectra, from Ref. [15] (first order in 
Hubble and sound flow parameters and first order in the UA) and Ref. [13] (second order 
in Hubble and sound flow parameters and first order in the UA), to Ref. [16] (second order 
in Hubble and sound flow parameters and third order in the UA), tells us that this is not 
an easy task. The enormity of the calculations involved, the algorithms to be designed to 
make it practically doable, makes such a task too complex to be included here; it deserves a 
separate study of its own®. 

The all-order finite quantities. We deal in this study with seven quantities namely, 
power spectra (2), tensor-to-scalar ratio (1), spectral indices (2), and running of the spectral 
indices (2). The first three exhibit logarithmic divergences while the last four do not. The 
working equation for the power spectra shares the same form. It follows that, should the 
logarithmic divergence appear in Pg, we expect it to also appear in Pt unless the divergence 
is only tied to the speed of sound which is clearly not the case here. The tensor-to-scalar 
ratio being the quotient of two divergent quantities, is likely to inherit such divergence; this, 
we find to be the case. Now, as can be traced in Subsec. 5.2, in comparison to that of 
Subsec. 2.3 for the power spectra, the remaining four quantities in our list share a different 
structure of the working equation making them compatible with the expansion scheme. In 
particular, in accord with the calculation in Sec. 5.2, with the placement of i^/9t) as 
the denominator of the integrand for Ug {nt), the end result is finite. Thus, these all-order 
hnite quantities (spectral indices and running of the spectral indices), apart from the possible 
imprecision associated with the preliminary exclusion of some correction factors in the UA, 
are valid. 

Origin of the divergences. The logarithmic divergences come from the incomplete can¬ 
cellation of logarithmic terms brought about by the expansion scheme (laid down in [13] 
and elaborated in Subsection 3.2) used in the calculation of some physical quantities such 
as the power spectra and tensor-to-scalar ratio, through the UA. Such divergences cannot 
be attributed to the prescription involving taking the limit^® as r/ —)• 0“ in the working 
equations (see for example, (2.21) and (2.23)) for the power spectra among others, because 
these working equations are mathematically well-defined in the mentioned limit. As such, 
the divergences are not a consequence of the UA itself but a mathematical artifact of the ex¬ 
pansion scheme used in the UA for the computation of the power spectra and tensor-to-scalar 
ratio. 

This implies that the expansion scheme has a characteristic that makes it incompatible 
with the working equation for the power spectra (and for the tensor-to-scalar ratio being just 
that—the ratio of the two power spectra). Such a characteristic cannot be attributed to the 
way we expand based on the variable of differentiation. Consider for instance, (3.14) for the 
first slow-roll parameter. 

e\ — ei* — ei*e 2 * in-h • • • 

rj* 

®There is a way (although also tedious,) to possibly disprove the lower order finite results for the power 
spectra : Compute dlnPs/dlnfc (or dlnPt/dlnfc) and compare the result with that calculated through the 
working equation for the corresponding spectral index. Note that although it can disprove the former if the 
two results do not match, it might not be able to prove the validity of the lower-order finite results for the 
power spectra. 

^'^More technically, it is kg 0“. However, k cancels out for every pair (fc, J?) or {k,T)\ for instance, in 
(4.7), Inw = ln^ = In^. 
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This end result is the same whether we use (3.11) involving N or the more computationally 
cumbersome (see Sec. 4) expansion exemplified by (4.25) involving rj: 


Cl — + 


dei 

dN 


(AiV) + 


1 d^ei 


- _L 

Cl — Cl* + — 

dr] 


iv - V*) + 


2!diV2 

1 d^ei 


2!diV2 


{ANf + ■■■ 
iv - V*f + 


However, note that many (in fact, infinite) combinations of the Hubble and sound flow 
parameters are of the same order; e.g., (ef^, ei*e 2 *, 62 * 63 *, • • •) ~ 0{e^). This gives us a hint 
that although the expansion for a given background quantity with respect to the Hubble and 
sound flow parameters does not depend on the variable of differentiation in the corresponding 
Taylor series, the expansion itself may not be unique. And if it is not unique and at the same 
time incompatible with the working equation for the power spectra, then the Hubble and 
sound flow parameters may not be the appropriate perturbation expansion parameter for the 
power spectra {Ps and Pt). The characteristic we are looking for, then, is the base of the 
expansion itself—the Hubble and sound flow parameters. 

We recall that unlike that of a typical perturbed equation of motion (e.g., that of a simple 
harmonic oscillator whose potential has an added small quartic term involving displacement 
of the oscillator about a certain equilibrium position), the Mukhanov-Sasaki equation does 
not contain an explicit perturbation parameter. As such, one should be careful about using 
any expansion parameter for any quantity whether it be physical or not, defined or derived 
from the Mukhanov-Sasaki equation. The suitability of an expansion scheme with respect to 
a given set of expansion parameters rests on the structure of a given (approximating) equation 
for a given quantity (in relation to the Mukhanov-Sasaki equation). For the case in hand, as 
already pointed out above, the current expansion scheme is incompatible with the working 
equations for the power spectra and compatible with that of the spectral indices. In relation 
to this, we note that there may be another group of perturbation expansion parameters with 
respect to which the power spectra is finite to all orders. Whether such finite results would 
agree and to what extent, with the finite lower-order results for the power spectra given in 
Sec. 4, remains to be seen. 

Physicality of the divergences. There are a multitude of ways by which a physical 
interpretation can be attached to the logarithmic divergences involved in this work. The 
expansion for the first slow-roll parameter for instance, 61 = 61 * — 61 * 62 * ln{r]/r]^,) + - ■ ■ , may be 
seen as an improper characterisation of the background. That this does not work for the power 
spectrum and works for the spectral index, seems to contradict this interpretation. We take 
the view that they are primarily mathematical in nature. Prematurely attaching a physical 
meaning may only complicate the process of solving the problem (in future work) brought 
about by the infinite final results. It is the mathematical aspect of the implementation of 
the UA that exhibits the problem and not the physics that it is trying to describe. 

Resolution in sight. Since the origin is mainly of mathematical nature namely, an artifact 
of an expansion scheme, the most probable resolution in sight should arise with respect to 
the same nature—find another expansion scheme. Such an expansion scheme should avoid 
as much as is (computationally) practicable, terms^^ that can potentially lead to a divergent 

^^Needless to say, we are not closing the possibility that a similar expansion scheme involving terms such 
as ln(77/?7*) may be able to lead to a finite result. In such a case, a complete cancellation of all potentially 
divergent terms is all that is needed. 


- 26 









final result. Moreover, it should retain the pragmatic decomposition of the background 
quantities (e.g., slow-roll parameter, scale factor, speed of sound, etc) in ascending order 
with respect to some set of parameters; see for instance, the expansion of ei given by (3.14). 
This should be contrasted with the (impractical) alternative expansion scheme exemplified 
by (4.25) where the higher-order derivative terms involve a mix of lower- and higher-order 
sub-terms in the Hubble and sound flow parameters. 


7 Concluding Remarks 

The uniform approximation (UA) is one of the leading semi-analytical approaches used to 
solve the Mukhanov-Sasaki equation and hence, to calculate relevant physical quantities in 
inflationary cosmology. It is a systematically improvable scheme with finite error bounds and 
does not require constant slow-roll parameters [22]. Nonetheless, the practical implementa¬ 
tion of this method at least for slow-roll ^-inflation, relies on a suitable expansion scheme in 
order to express physical quantities in terms of some possibly measurable parameters. 

We note the expansion scheme for slow-roll parameters and other background quantities 
(e.g., scale factor, speed of sound, etc) introduced in [15], expounded in Ref. [13], and adopted 
in Refs. [14, 16]. This expansion scheme involves powers of ln(? 7 /r/*). With the prescription 
that the limit as r/ —>■ 0 “ has to be evaluated at the end of the calculation, this may potentially 
lead to logarithmically divergent result. In their calculation of the power spectra up to second 
order with respect to the Hubble and sound flow parameters, the authors of Ref. [13] stated 
a “consistency check” as to the validity of their implementation of the UA that can be stated 
as : all terms that can potentially give rise to logarithmically divergent result should exactly 
cancel prior to taking the limit as 77 —?• 0“. Their second-order calculation successfully 
passed this “consistency check,” but the all-order check was not carried out. There was also 
no rigorously sufficient characterisation of all forms of divergences that may come out of 
the calculation. (This seems unnecessary as long as there is an assurance that all forms of 
divergences will cancel out in the end; it turns out not to be the case.) 

This work is an attempt to fill in these gaps. Focusing on slow-roll /c-inflation, we (a) 
characterised the form of the divergences as purely logarithmic in nature and (b) set the limit 
as to how far one can calculate several physical quantities in terms of the Hubble and sound 
flow parameters before hitting a logarithmically divergent result. It turns out that the power 
spectra for both the scalar and tensor perturbations can only be calculated up to second 
order in terms of the Hubble and sound flow parameters (see Table 1). The tensor-to-scalar 
ratio is likewise bounded up to third order. However, the spectral indices and the running of 
the spectral indices can be calculated to any order. So far, the most precise calculations of 
these physical quantities in fe-inflation through the UA can be found in Ref. [16] where the 
calculation is carried out up to third order of precision in the UA. The boundary that we set 
here is consistent with this work. 

Looking back, we find nothing inherently wrong with the mathematically well-defined 
working equations provided by the UA for the power spectra and tensor-to-scalar ratio. It is 
with the practical implementation of these working equations that we get the logarithmically 
divergent results in the higher-order parts. Such divergences stem from the use of the expan¬ 
sion scheme involving potentially divergent terms that do not find a complete cancellation 
in the end. It turns out that with respect to the structure of the working equation for the 
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Physical Quantity 


0 12 3 4 


power spectrum (scalar), Pg o o 
power spectrum (tensor), o o 


tensor-to-scalar ratio, r o 

spectral index (scalar), Ug o o 

spectral index (tensor), n* o o 

running (scalar), ag o o 

running (tensor), at o o 


o • 
o • 

o o 
o o 
o o 
o o 
o o 


o 

o 


Table 1. Logarithmic divergences in the calculation of the power spectrum and other physical quan¬ 
tities. Filled circles represent results involving logarithmically divergent terms while unfilled circles 
represent finite results. The numbers in the first row indicate the order of a given physical quantity 
with respect to the Hubble and sound flow parameters. Note that once a logarithmic divergence 
appears (as in the power spectra and tensor-to-scalar ratio), it permeates to all succeeding orders. 
The spectral indices and running are finite to all orders. 


seven quantities listed in Table 1, the expansion scheme is incompatible with that of the 
power spectra and tensor-to-scalar ratio and compatible with that of the spectral indices and 
running of these indices. The calculation results for the last four quantities {ug, nt, Og, at) 
are therefore, completely valid apart from possible imprecision that can be remedied by in¬ 
creasing the level of precision of calculation in the UA. The first three quantities on the other 
hand, do not rest on a solid mathematical foundation. For this reason, some of the terms 
in the lower-order finite results for Pg, Pt, and r may be misleading. However, owing to the 
enormity of the calculations involved needed to verify this, we set aside the necessary analysis 
of these lower order terms for future studies. 

We hope that we have clarified the limits of one significant implementation of the UA 
in the calculation of several important physical quantities in inflationary cosmology. Ours 
is a rigorous identification and characterisation of the problem in relation to the existence 
of logarithmic divergences in an otherwise, finite set of quantities. In our view this puts the 
current works in the literature into a proper perspective. Moreover, this sets up the framework 
of how to solve a problem that has only just been formally identified and characterised in 
this work: How do we implement the UA in the calculation of power spectra for slow-roll 
A:-inflation without suffering from logarithmic divergences to all orders with respect to a given 
set of perturbation expansion parameters? 
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